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1. Introduction 

The recent experimental evidence of the quark-gluon plasma phase of QCD in RHIC pro- 
vides experimentalists and theorists with a formidable challenge: How is the physics of this new 
phase uncovered from the products of the fireball when these products have necessarily condensed 
back into the confined phase? 

To answer this question it is crucial to know whether hadronic states in fact persist in the 
quark-gluon plasma phase. Current understanding of the plasma phase suggests that while it is not 
confining, it is still strongly interacting, suggesting that hadronic states can indeed survive above 
the de-confining temperature, T c [Q]. 

This study investigates mesonic states at and above T c via their spectral functions, p(ft),p). p 
can be defined from Euclidean hadronic 2-point function, G(?,p), as 

r°° dco 
G(t,p)= J p(a,p)K(t,(o)— (1.1) 

where the kernel, K(t, co), at temperature, T, is given by 

cosh[(»(f-l/(2r))] 



K{t,co) 



sinh[fl)/(2r)] 



Spectral functions contain information on the stability of hadronic states, transport coefficients 
and dilepton production rates, and thus are fundamental physical quantities of the system. They 
can be derived from lattice calculations of hadronic correlators via the Maximum Entropy Method 
(MEM)[g]. However, the success of this approach relies on the quality of the input data, in partic- 
ular on having estimates of G(t) over a large time range. Using isotropic lattices, this requirement 
runs at odds with finite temperature studies which restrict the temporal extent to be small. We 
overcome this conflict by using anisotropic lattices, to ensure that the number of time points, N t , 
is large, while maintaining a non-zero temperature, T = \/{a t N t ) (where a t is the temporal lattice 
spacing which we take much smaller than the spatial spacing, a s ). Using this approach, we have 
performed dynamical studies of charmonium states at temperatures T c <T < 2T C and are able to 
confirm that they are still bound even for temperatures above T c . This work is a continuation of that 
presented in [Eh. 



2. Lattice Parameters 



Because our lattices used dynamical fermions, tuning the two anisotropy parameters in the 
gauge and quark actions is highly non-trivial. This is because there is feedback from the fermions 
to the gluons which is not present in the quenched anisotropic case. Details of this procedure for 
the choice of light quark mass and gauge coupling used in this work can be found in [Q]. The 
ensembles studied are more highly tuned than in our previous work presented in [Q] and have a 
renormalised anisotropy for both quarks and gluons very close to 6. 

A highly improved anisotropic gauge action was used with Wilson + Hamber-Wu fermions 
with stout links. The details of this action can be found in [Q]. Table |l| lists the parameters used in 
the simulations. As can be seen, three temperatures were studied: T/T c ~ 1, 1.3, 2 and two spatial 
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Light quarks 


Mn/Mp 


~0.5 




Heavy quarks 


am c 


0.080,0.092 




(Renormalised) Anisotropy 




6 




Lattice spacings 


a, 


~ 0.025 fm 






a s 


~0.15fm 




Spatial Volume 


N 3 


8 3 (&12 3 ) 




Temporal Extent 


N 


16 < — ► 


~ 2T C 






24 < — ► 


~ 1.3 T c 






32 < — ► 




Statistics 


Nfg 


~500 





Table 1: The lattice parameters used. 



volumes were simulated enabling a finite volume analysis. Two heavy quark masses were studied, 
with the am c = 0.080 being closer to the experimental charm quark. 

The deconfining transition was found between N t = 34 and 33 by studying the Polyakov line 
on the Nf = 12 3 spatial volume. (No discontinuity in the Polyakov line was found in the Nf = 8 3 
case.) 

Before embarking on the full simulation, it is worth studying the systematic effects induced 
by the lattice discretisation. Figure [j] shows the pseudoscalar and vector spectral functions for the 
free theory for both the continuum and lattice. (See ^ for details of how this calculation was 
performed.) In the lattice case, the anisotropy of 6 was chosen and two temporal extents, N t = 24 
and 32 are shown. As can be seen, the lattice functions reproduce the continuum ones for small 
and moderate energies, but there is an unphysical cusp at a t co ~ 0.7 which is independent of N t and 
channel. This leads us to conclude that lattice spectral functions above this a t CO value are likely to 
be contaminated by systematic effects. Note also that the cusp for the run parameters in [Q] occurs 
at the smaller value of a t co ~ 0.6, and hence the results presented here are closer to the continuum 
ones over a slightly larger energy range. 

3. Results 

Since we have correlation functions calculated for a number of different parameters, we can 
study how the spectral functions vary with: 

• energy resolution for the MEM 

• starting time for the MEM analysis 

• spatial volume 

• heavy quark mass 

• channel (i.e. pseudoscalar, vector, axial, scalar) 

• temperature 
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fl ? CO 

Figure 1: The spectral function for the vector and pseudoscalar mesons in free field theory for the parameter 
values in this work. The continuum results and those from the lattice with N t = 24 and 32 are shown. 




Figure 2: LEFT: The spectral function with various energy resolutions, i.e. the number of discretised 
CO points, N a , showing the insensitivity of the spectral function to this quantity. The scalar channel with 
am, = 0.092 on the 8 3 x 24 lattice was used. RIGHT: The spectral function with various starting times, 
r start f° r tne MEM analysis. This shows the insensitivity of the spectral function to f s tart within the range 
2 < fstart < 4. (The same channel as the left panel was used.) 

The hope is that the spectral function determined from MEM will not vary with the first three 
items since they are unphysical, and this is what we find. The left panel in figure |2| shows the total 
lack of effect when the energy resolution of the MEM is varied from n a = 100 through to 1000. 

In the right panel of figure |2| the effect of varying the start time, f s tart> of the MEM is studied. 
Note that the spectral functions for ? s tart = 2,3 or 4 are all compatible, i.e. there is a range of 
sensible £ s tart values which produce consistent results. However, for £ s tart = 1 an d ^start > 5 the 
spectral functions do vary. This can be understood since there are contact terms which enter the 
? start = 1 case, and if f s tart > 5 is used, there is little signal left for the MEM to latch onto (N t = 24 
in this case). 
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Pseudoscalar 8~ x24 Pseudoscalar 




CO 



Figure 3: LEFT: The spectral function for both spatial volumes N s = 8 and 12 illustrating the small finite 
volume effects. The pseudoscalar channel with N t = 16 is shown. RIGHT: The spectral function with both 
am c = 0.08 and 0.092 showing the small sensitivity to the m c value. The pseudoscalar channel on the 8 3 x 24 
lattice was used. 

Finite volume effects are studied in left panel of figure ||. As can be seen, there is little variation 
in the spectral function as N s is changed, confirming that our results are not dominated by finite 
volume systematics. 

Now that these lattice systematic effects are studied and have been shown not to contaminate 
our spectral functions, we turn to varying the physical parameters. The right panel of figure || 
shows the effect of varying the charm quark mass from am c = 0.08 to 0.092. The ground state 
peak is slightly higher in the latter case as expected, but otherwise the spectral function shows little 
dependency on m c . 

Finally we consider varying the temperature. Figures ^| show the spectral functions for the 
three temperatures studied for the four mesonic channels. (N s = 8 is used throughout.) As can be 
seen there is clear evidence of a bound state in the pseudoscalar and vector channels which survives 
above T c and melts between 1.3 T c and T c . At 2T C there is evidence of non-zero spectral weight at 
small energies. In the case of the axial channel, the state seems to melt by around ~ 1.3 T c , and for 
the scalar channel the situation is less clear. 

4. Conclusions 

This work continues that first presented in [|J] where anisotropic dynamical lattices were used 
together with the MEM technique to determine charmonium spectral functions at non-zero temper- 
atures. Improvements over the work described in [Q] were that the anisotropy parameters are now 
better tuned, a larger volume considered with more than one heavy quark mass simulated. 

We have confirmed the spectral functions from MEM to check that they are stable for sen- 
sible variations in the MEM parameters (energy discretisation, starting value of the time interval, 
volume). 

Our results confirm our earlier work and previous quenched studies [^] : the pseudoscalar 
(T7 C ) and vector (J/y) charmonium states melt between 1.3T C and 2T C , while the scalar (Xm) an d 
axialvector (% c i) states are less less certain, but seem to melt by 1.3 T c . 
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Pseudoscalar Vector 




Figure 4: The meson spectral function at various temperatures (with am c = 0.08 and N s = 8). 

Future work will include varying the default model, inclusion of non-zero momenta and studies 
of light hadrons. 
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